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The widths and overall profiles of dielectric grating lines can be determined by measuring the 
intensity of diffracted laser light from the sample over a specified range of incident beam angles. 
This techruque. known as 2-9 scatierometry, is able to accurately and precisely measure photoresist 
structures in the subhalf micron regime; Moreover, a 2-0 scatterometer is capable of making 
measurerhents iii a rapid and nondestructive manner. To test this technique we measiired five , 
identically processed wafers with nominal 0.5 /un line/0.5 /xm space grating patterns. Each wafer 
comprised gratings in a Shipley 89131 negative photoresist exposed in a matrix of incremental 
exposure doses and focus settings. The scatterometry results were consistent with cross-sectional 
and top-Kiowri scanning electron microscopy (SEM) measurements of the sanie structures. The i 
average deviation of 1 1 scalteroineter linewidth measurements from top-down SEM measarements, 
over a broad exposure range, is 14.5 nm. In addition, the repeatability (la) of the 2-0 scatterometer 
is shown to be excellent: 0.5 nm for consecutive measurements arid 0.8 nm for day-to-day 
nieasurements. © 1995 American Vacuum Society. 



I. MICROELECTRONlbS METROLOGY: AN 
OVERVIEW 

As the ihtegratioii and speed of microelectronic devices 
increase, the dimensions of circuit structures must continue 
to decrease in size and to improve in terms of profile edge 
acuity. High-end devices such as microprocessors require a 
considerable number of process steps. Therefore, it is becom- 
ing increasingly important to have an accurate, quantitative 
description of the submicron structure present at each step. 
This places incr^ing demands on present-day metrology 
techniques, and, at present^ there is no single technique ca- 
pable of meeting all of these demands. Fuithermore, there is 
'a growing need for wafer process monitoring and closed loop 
control, both of which would require in-line sensors. In ad- 
dition, these wafer/process sensors must bie integrated into 
the fabrication line, i.e., they must have a reasonable foot- 
print and be able to monitor the process rapidly. Most line- 
width measurement equipment is either too slow or too com- 
plex to be implemented in an in-line arrangement for high- 
throughput wafer evaluation. An ideal in-line metrology 
technique or sensor would be capable of making measure- 
ments that are sufficiently accurate, repeatable, and rapid at a 
low cost. 

The major tools at the disposal of the semiconductor in- 
dustry include electrical measurements, optical microscopy, 
various probe microscopy techniques, and scanning electron 
microscopy (SEM). Each technique has its strengths and 
weaknesses, and many in industry believe that these metrolo- 
gies will be used in a complementary fashion to exploit the 
potential of each. 

Electrical testing is a useful method for testing submicron 
features under certain conditions. The most appealing at- 
tributes of electrical methods include 6 nm precision (3a)^ 
and the ability to measure critical dimensions as small as 0.2 
/on. A major limitation is the fact that measurements must be 



made oh a conducting filin after the etching step when re- 
working the wafer is less feasible.''^ In addition, because it 
involves contactj it is slow and not araehablc to in-line ap- 
plications. " . ; ■ 

Optical microscdpiy still plays an important role in^ubini- 
cron metrology due to its versatility, simplicity, and speed.* 
Sectioning optical thicroscopes have received much attention 
lately for achieving three-dimensional imaging of siibmicfon 
structures, an attractive feature for overlay measurernent.^ 
The two predominant types of sectioning optical micro^ 
scopes are confocal sbanning optical microscopes (CSOMs) 
and interference effect microscojpes.^ One of the most attrac- 
tive features of optical microscopes is that they do not re- 
quire placing the sample in a vacuum, a fact that translates 
into lower cost relative to niost competing nietrologies. Uii- 
fortunately, traditional microscope designs suffer in terms of 
resolution and precision when measuring structures at or be- 
low 0.5 /xm.** Improvements in technique, design, and inate- 
rials, however, may continue to improve accuracy and re- 
peatability and thus keep optical microscopy a viable 
metrology. Diffraction limitations are driving optical micro- 
scope design into scanning near-field optical microscopy 
(SNOM)'*'^ which allows for high-resolution imaging. It is 
important to note that SNOMs can also be classified as probe 
microscopes,* another metrology for making three- 
dimensional measurements. 

Although in general probe-type instruments are capable of 
high resolution in three dimensions, the measurement pro- 
cess is slow. Furthermore, the image obtained is difficult to 
interpret unless the tip shape is well known.* Probe tips miist 
be highly specialized for the type of topography to be 
scanned and, in some cases, especially when high aspect ra- 
tio features are involved, portions of the profile can be inac- 
cessible to the probe. Scanning tunneling microscopes 
(STMs) are a special class of probe microscopes capable of 
atomic resolution in a high vacuum. There are two major 



1484 J. Vac. Sd. TechnoL B 13(4), Jul/Aug 1995 0734-211 X/95/13(4)/1 484/1 2/$6.00 ©1995 Arribrican Vacuum Society 1484 



1485 



Raymond ef a/.: Metrology of subwavelength photoresist gratings 



1485 



drawbacks to direct-current tunneling nniicroscopes; they can:: 
not measure nonconductive samples and close contact be- 
tween the probe and the sample results in forces that can . 

waiter both.*;.; . 

v i. Samples measured with atorriic force- microsc6j)es. 

; (AFMs) need not be coated with a metal because no current 

' flows from the probe to the sample. Atomic resolution has 
been reported for AFMs operating in repulsive mode,^'*^ 

■■ where the probe-sample distance is on the order of I mn or 
less. As the probe-sample separation increases,^ probe- 
sample forces become attractive and resohitidn generally de-. 

, creases. Recent literature reports lateral resolutions in the 
5-10 hm range.'* One drawback to using an AFM is that the 
probe tip must be frequently and accurately characterized, a 
tedious process that is not always practical. 

SEMs, which offer greater resolution than optical micror 

■ scopes, are becoming the standard metrology technique for 
the semiconductor industry.'^ The enhanced resolution 
comes at a price, however; sarnples must be measured in a 
high vacuum, thereby increasing the complexity and cost of 
each measurement and limiting throughput.* In the case of 
samples with insulating structuresi charging effects and elec- 
tron absorption can cause measurement eaors and sample 
damage^ respectively. 

; Top-down^ or line-scan (LS) SEM measurements can pro- 
vide linewidth infdrnlation with the use of an appropriate 
edge detection algorithm. Although the use of different edge 
detection algorithms has lead to inconsistent results," recent 
improvements in electron-sample interaction models are en- 
abling researchers to make more accurate characterizations. 
By incorporating backscattered electron (BSE) models and 
Monte Carlo algorithms, edge positions have been located to 
within 10 nm, and sensitivity to rough and rounded profiles 
has bfeen demonstrated. Backscattered electron imaging has 
also been reported capable of measuring high aspect ratio 
resist features with a 3a- precision of 4 nm.*^ The calibration 
of Sl^s for linewidth measurement, however, has been dif- 
ficult due to the lack of availability of a standard.*^ Further- 
more, top-down linewidth measurements on nonsquare pro- 
files that do not have a well-defined vertical edge, can be 
misleading.** 

Scanning electron microscopes can also be operated in 
"cross-section" (CS) mode. While this approach allows for 
complete line shape characterization, studies have indicated 
that the systematic error of a CS scan is larger than diat of a 
top-down SEM measurement under the same operating con- 
ditions. This is thought to be a consequence of different con- 
trast formation mechanisms that result from the different 
sample geometries of the LS and CS modes. The distortions 
introduced in CS mode can cause errore that transcend the 
linewidth tolerances permitted for subhalf micron 
structures.*' Furthermore, cross sectioning obviously de- 
stroys die wafer, a drawback that rules this technique out for 
high volume, in-line applications. 

I>iffraction-based analysis techniques** such as scaiterom- 
etry arc especially well suited for microelectronics metrology 
applications because they are nondestructive and rapid. *^ 
These metrology techniques have already proven useful in a 
variety of microlithographic applications, including 
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ahgnment,^ overiay,^ • temperature measurcraent,^ monitor- 
ing of latent image focus,^ exposure,^'^ and postexposure 
bake.^ Scatterometry has shown promise as an accurate, re- 
peatable, aiid inexpensive alternative to other metrologies.*^ 
With its sensitivity, versatilityi and capability of high- 
throughput measurement; scatterometry is gaining recogni- 
tion as a potentially powerful metrology of the very near 
future. 

IL SCATTEROMETRY 

Scatteronietry is the angle-resolved nieiasuremeii^ and 
characterization of light scattered from a structure. For struc- 
nires that are periodic, incident light is scattered or diffracted 
into different orders. The angular location ^ of the mth dif- 
fraction order with respect to the iangle of incidence 6^ is 
specified by the grating equation 

sin ^,+sin Sf=m—\ (i) 

where X is the wavelength df incident light and d the period 
of the diffracting structure. 

The diffracted light pattern from som6 structure can be 
used as a "fingerprint" or "signamre** for identifying the 
dimensions of the structure itself. In addition to period, 
which can bei determined quite easily, more specific dimen- 
sions such as the thickness of thd photoresist, angle of the 
sidewalls, arid width of the line can also be measured by 
analyzing the iscatter pattern. 

The scatterometric analysis can best be defined in two 
steps. First, the diffracted light fingerprint is measured. This 
is known as the forward problem. The second step, the in- 
verse problem, consists of detenriining the shape of the grat- 
ing that scatters the incident light. To solve this problem, the 
grating shape is parameterized^ and a parameter space is 
defined by allowing each grating shape parameter to vary 
over a certain range. A rigorous diffraction modeP^ is used to 
calculate the diffracted light fingerprint from each grating in 
the parameter space, and a statistical prediction algorithm is 
trained on tiiis theoretical calibration data. Subsequently this 
prediction algorithm is used to determme the parameters that 
correspond to the fingerprint found in the forward problem. 

In our work discussed here, the forward problem is solved 
using a 2-0 scatterometer. seen in Fig. 1. A simple He-Ne 
laser beam, X=633 nm, is incident upon a sample after it 
passes through a spatial filter and some focusing optics. The 
sample is mounted on a stage that permits it to rotate. Since 
the beam itself is fixed, this rotation changes the angle of 
incidence on the sample. Using grating Eq. (1), the detector 
arm of the scatterometer is able to follow any diffracted or- 
der as this angle is varied. The intensity of a particular dif- 
fraction order as a function of incident angle, which is 
known as a 2-G plot or scattier "signature," is then down- 
loaded to a computer, where the analysis can begin. 

Previous scatterometry applications on large period 
samples measured the entire diffraction order envelope at a 
fixed angle and used this as a means for determining the 
dimensions of the scattering structure.^*^ But, for I /un 
period gratings, die size of interest for this study, only the 0th 
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and ±lst diffraction orders exist over a practical angular 
range. Measuring this dififraction envelope, then, does not 
provide the data required for an accurate analysis. However, 
a 2-9 plot contains enough information for proper character- 
ization of the diffracting structure and may be thought of as 
a fingerprint of that structure. For example. Fig. 2 shows two 
2-9 plots of the zeroth order intensity in which the linewidth 
of a 1 /im pitch grating is changed by 10 nm. 

Hie scatterometer*s laser was TE polarized for this study. 
To enhance the signal-to-noise ratio and minimize error, a 
conventional lock-in detection system is used. The beam spot 
size on our sample was approximately 400 /xm for ease of 
use, but can be made smaller, and if changed, does not affect 
the diffraction efficiency.^"^*^^ One important consideration 
regarding spot size is the number of grating periods illumi 
nated by the beam. Theory assumes the structure is infinitely 
wide, in practice, 20 or so illuminated period are sufficient to 



approximate an infinite structure.^** Thus, for a 1 pirn pitch 
grating, the beam must have a diameter of at least 20 ^un. 

It is important to note a major difference between meia- 
surements made by scatterometry and those made by other 
methods. Because the spot size of the incident beam is much 
larger than die period of the diffracting structure^ many lines 
and spaces are illuminated and all contribute collectively to 
the diffraction signal. Hierefore, line shape characterization 
values using scatterometry are averages of all the illuminated 
structures. This averaging effect can be an advantage because 
a SEM photo may depict a structure that is not representative 
of the entire sample. The spot size of the scatterometer can 
be adjusted to measure the average line shapes of hundieds 
of features or a localized group of therh. This illustrates how 
competing technologies are often complementary in certain 
respects. 

fll. DiFFRACmbri MODEL 

With the acquisition of experimental diffraction data us- 
ing the 2-0 scatterometer, the forward problem is complete. 
The next step in the overall analysis strategy is to generate a 
series of dieoretical diffraction data so that the inverse prob- 
lem may be solved. This requires a diffraction model.^^ 

For structures which arc periodic in one dimension, these 
theoretical diffraction "signatures" can be constructed using 
the rigorous coupled wave tiieory (RCWT) developed by 
Gaylord and Moharram.^^ This vector differential technique 
uses Maxwell's equations and the appropriate boundary con- 
ditions at all interfaces to solve for die electric fields in each 
region of the sample. The total reflected field intensity can 
then be calculated; it is comprised of discrete diffracted "or- 
ders" as a result of the periodic structure, as is evidenced by 
grating Eq. (1) presented earlier. The RCWT diffraction 
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; model is particularly useful for photoresist grating applica- 
tions because it can account for nonsquare grating profiles, 
discussed later in Sec. VI, as well as any number of under- : 
lying film layers or substrates. A complete overview of rig- 
orous coupled wave theory with respect to scatteroraetric ap- 
plications can be found in Naqvi et al?^ 
■ The generation of theoretical scatter signatures is the most 
; time-consuming process of the entire scatterometric analysis. 

For example, it takes approximately 30 s for a Hewleft- 
. Packard (HP) Apollo workstation to generate one 2-B signa- 
; . ture for a square grating profile. Most of the analyses pre- 
sented here required some 400 signatures, resulting in 3.33 h 
of computation time. It is important to note that, for a given 
parameter space, the calibration procedure only needs to be 
performed oncc; 

IV. STATISTICAL ANALYSES/PREDICT^ 
ALGORfTHMS 

Once the theoretical database hais been generated, the sec- 
ond part of the inverse problem, namely^ characterizing the 
diffracting structure, begins. The partial least squares (PLS) 
algorithm of Chemometrics, a multivariate analysis package, 
was selected as one means for solving the inverse problem 
based on its success in the past^^'^^ and its good performance 
over a wide range of experimental conditions.^^ A detailed 
description of Chemometrics is given by Haaland et dl^^ 

With PLS, the statistical solution of the inverse problem 
occurs in two phases. In the calibration phase, either empiri-- 
cal or theoretically generated diffraction order intensity data 
are used to train a prediction model. All of the results pre- 
sented in this article involve the prediction of line shapes; 
and thus it was possible to generate the theoretical data sets 
using RCWT. If we had wished to predict stepper settings 
such as exposure dose 6r focus setting, measurements on 
gratings exposed at each value in the focus-exposure matrix 
woulSJiave been taken. These data Would then be used as the 
calibration set to predict focus and exposure dose on un- 
known samples. In the second phase, the parameters of in- 
terest of the unknown sample are predicted using the model 
derived from the calibration set. The prediction step is very 
rapid, allowing for a nearly instantaneous characterization of 
the unknown sample. ■ 

Another technique for solving the inversie problem is 
through the use of minimum mean square error (MMS) 
analysis. Like Chemometrics, the foundation for MMS 
analysis is in the calibration set, which, again, can utilize 
either empirical or theoretical data for training. The actual 
analysis proceeds by comparing an unknown experimental 
intensity data file with each file in the calibration set. The 
mean square error is calculated for each comparison, and the 
closest match (minimum mean square error) recorded. The 
linewidth or overall line shape associated with the closest 
calibration file is taken to be the parameters of the unknown 
structure. 

Both PLS and MMS have performed well in our analyses. 
The method of choice depends on the application involved. 
Typically MMS works better when intensity data are only 
weakly affected by a change in a given parameter. One ex- 
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ample is the slight Change in diffracted power caused by 
slightly rounded comers on photoresist lines. In addition, the. 
actual minimum mean square error returned by the algorithm . 
is a direct indicator of how well theory and experiment agree 
and this error value serves as a confidence level of the pre-: 
diction. However, MMS is riot capable of interpolating be- 
tween calibration data points as is PLS. Still, MMS has 
proven to be more dependable if the computing power is 
available to create a calibration set oif sufficient resolution. 
Another potential pitfall in using MMS is that it is not suf- 
ficiently sophisticated to recognize data outliers; it simply 
determines the difference between theory and experiment at 
each data point and does not account for similarities and/or 
differences between twa sets of data. S : 

Chemometrics/PLS actually builds a calibration model 
based upon an input calibration data set. Using this calibra- 
tion model, it is able to notice trends in data and predict an 
unknown parameter, even to the point of interpolating be^ 
tween calibration data based on those trends. Chemometrics 
also has built-in statistical returns that serve as confidence 
indicators for the prediction. However, when parameter 
changes manifest themselves only slighfly in the scatter sig-; 
nature, Chemometrics has not performed well. 

We have also evaluated the impact of the range of inci- 
dent angles on the effectiveness of the; prediction. In prin- 
ciple the 2-6 scatterometer can take data on any diffraction 
order for incident angles ranging from —90 to +90°. In prac- 
tice, however, we have limited the angular range to ±60°; at 
angles much beyond d[iis the incident light is almost com- 
pletely reflected and contains little information about the 
structure (i.e., the zeroth order intensity approaches simple 
Fresnel reflectances for near-grazing incidence). In the re-, 
suits presented here two angular ranges are lised: one from 
T- to 60^, the other from 2° to 42° (the 0^ and T angles 
cannot be measured because the detector occludes the inci- 
dent beam in the current experimental arrangement). The 
smaller range was studied because it corresponds to a range 
that could be covered by a prototype conmiercial scatterom- 
eter currently under development. 

The experinients discussed here conceiitrated on the use 
of zeroth order data for the analyses, although the use of the 
± 1st orders was explored as well. As will be seen, the zeroth 
order signatures alone contain sufficient information about 
the unknown structure to yield accurate characterization oif 
the diffracting structure. From an experiraetital standpoint 
using a single order for prediction purposes is advantageous 
because it minimizes the number of measurements. In addi- 
tion, as the microelectronics industry pushes toward smaller 
device dimensions, the ± 1st orders will eventually no longer 
exist for He-Ne illumination, making the 0th order the only 
choice for scatterometric measurements using 633 nm illu- 
mination. 

V, DIELECTRIC GRATING UNEWIDTH 
MEASUREMENTS 

Determination of the width, 6ften referred to as the criti- 
cal dimension (CD), is the first step in the scatterometric 
characterization of a subnaicron line in developed photore- 
sist. If shape parameters odier than linewidth are subtle to the 
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point of hieing irrelevant, this technique can provide an ex- 
cellent measurement of the near-rectaiigular profile. Results 
for 1.5 lines in 32 /mi period etched silicon structures 
have been published by Naqvi et al}^ Large period gratings 
have the advantajge of producing a large number of diffrac- 
tion orders, whereas gratings with a 1 fim period will diffract 
light into only one or two orders. 

Because they contain a wide variety of diffracting struc- 
tures, gratings generated within a fociis/exposure (F/E) ma- 
trix provide a very good test for scatterometry ajpplicationis. 
Changes to exposure, for example, predominantly affect the 
liiiewidth while focus changes generally influence the overall 
shape of the resist. The linewidth measurement study was 
initiated with one pilot wafer with a F/E matrix printed on it; 
this wafer is referred to as "wafer 3." Altfiough bodi the 
height and index of the resist was known to within a particu- 
lar range, the degree of uncertainty for each was unaccept- 
able for our model; Thus it was necessary to iise the scatter^ 
ometer to determine the height and index of the photoresist 
before predictions of any line widths could be performed. 

These two parameters can be determined fh)m either the 
diffraction data or 2-0 scans of an unpattemed region of the 
photoresist (i.e., a thin film). For a 2-6 scan on a thin film, 
the reflected intensity goes through a minimum whenever 
destructive interference occurs. By generating a theoretical 
database of index/height combinations and comparing them 
via statistical algorithms to experimental scans of the photo- 
resist, we were able to determine the index and height of the 
resist layer. The results matched the specified nominal val- 
ues: the height was determined to be 1.00 /xm, and the re- 
fractive index 1.576. In addition, from examining a pat- 
terned, diffracting sample and using zerddi order diffraction 
data, the MMS algorithm has been able to do a three param- 
eter search for height, index, and linewidth: This technique 
determined the heigjit to be 0.990 fim and the index 1.575. 

With the two photoresist parameters determined it was 
f^^ssible to proceed with the analysis of grating linewidths on 
wafer 3. Initially we believed the sidewalls of the photoresist 
lines would not be perfectly square and would have some 
angle to them. But, when early predictions were run with the 
sidewall angle as a free parameter, we found the angles were 
being predicted at 90*", i.e., the developed lines had vertical 
sidewalls. Thus it was only necessary to keep one parameter 
free in the final prediction model: liflewidth. For wafer 3 the 
prediction model linewidth range ran from 0.350 to 0.750 
Mm in 1 nm steps. 

After the predictions were performed, comparisons be- 
tween cross-sectional SEM and scatterometry measurements 
for wafer 3 showed good agreement. Unfortunately the SEM 
measurements were limited to six locations, three of which 
were useful for linewidth measurements and 3 for profile 
determination. These profile results are discussed in a later 
section. But the real impact of this one wafer pilot study was 
that it launched a more thorough and complete study of the 
scatterometry technique, with five wafers and approximation 
60 top-down and 12 cross-sectional SEM measurements. 

To further test the scatterometer and linewidth prediction 
algorithms, five sample wafers, labeled 4-8, were used. 
Each wafer had grating structures of 1 /mi pitch and 0.5 /tm 
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nominal linewidth. These griatings had been exposed and de- 
veloped in photoresist on a silicon substrate for a range of 
focus arid exposure values. The exposure values extended 
from 30 to 100 nJ in 5 inJ increments, and the focus ranged 
from 0 to ±1:4 /^m in ±0.2 /tm steps. The 2-9 scans Wett 
performed on each of the wafers, with data collected froiii 
both the 0th and ± 1st orders. The data was then used in the 
prediction algorithms to determine the linewidth of die grat- 
ing. Hie calibration set , was generated using the RCWT 
model, with the assumptioiJ that the overall grating shap^ 
was rectangular. 

The effects of using different linewidth prediction algo- 
rithrtis for each of the wafers can be seen in Figs. 3-6. For 
comparison purpose, SEM measurements are included in 
some of the plots. Figures 3 and 4, which show the results 
for wafer 5, indicate very good agreement between both top- 
down and cross-sectional SEM and scatterometric measure- 
ments. The MMS algorithm: predictions (Fig. 3) are envel- 
oped by the top-dowri and cross-sectional SEM 
measurements over the entire exposure range. Tlie PLS algo- 
rithm results for this wafer can be seen in Fig. 4; the results 
are not as close as those using MMS (particularly at low 
exposure energy), but they are still consistent with the SEM 
measurements. With regard to angular range, the 60** and 42** 
datasets for either statistical algorithm yield approximately 
the same linewidth — the most significant difference between 
the two occurs at the 85 mJ exposure dose. 

Results for wafers 7 and 8 (Figs. 5 and 6) demonstrate 
less agreement with top-down SEM measureihents over the 
entire exposure range. However, for both wafers all of the 
results show a clear offset between the SEM and 
scatterbmetry-measured dimensions, regardless of the predic- 
tion algorithm used. For example, the standard deviation 
(la-) of the difference between SEM and scatterometry mea- 
surements (42° MMS algoridim) for wafer 7 is 7.9 nm. For 
wafer 8, the same algorithm results in a deviation of 8.2 nm. 
So, despite the offset, the scatterometry predictions are very 
consistent with the SEM measurements. Differences between 
the two techniques could perhaps be due to the effects of 
subtle features on the resist profiles. In other words, the theo^ 
retical model is assuming that the profiles are exactly square 
when, in fact, there are subtle features on the resist line that 
affect the scatter signatures. These complex features will be 
accounted for in Sec. VI. 

Finally, the results for wafer 4 can be seen in Fig. 7. For 
this wafer two sets of top-down SEM measurements were 
performed. The average differeince between the two SEM 
measurements is 12 nm. The results show approximately the 
same amount of offset as that for wafers 7 and 8 and are very 
linear with the SEM values. This can be seen in Fig. 8, which 
is a plot of the SEM measurements versus the scatterometry 
results. The correlation coefficient for the linear fit between 
the two techniques is 0.9922. Using the MMS 42° algorithm 
in comparison to SEM "b," the standard deviation (la) of 
the difference between the two measurement techniques is 
7.4 nm. The average differences and standard deviations 
from SiEM measurements for wafers 4, 5, 7,' and 8 are sum- 
marized in Table 1, 

The results presented thus far have all b^eh obtained us- 
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Fig. 7. The PLS and MMS scatterometry results for wafer 4, shown with two SEM measurcment«i on ihe same location. 
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Fig. 8. The. scatterometry results for wafer 4 vs the SEM measurements. 
JVST B - Microelectronics and Nanometer Strtictures 



0.8 



i49f2 Raymond et aU Metrology of subwaveleiigth photoresist gratings 



1492; 



Table I. Summary of average differences and standard deviations. Campari- 
scm of scatterometiy measurements (MMS 42**) to top-down SEM results; 







Standard deviation 


Wafer 


Average difference 


(1 a) of difference 


LD. No. 


from SEM values 


from SEM values 


4 


54.0 nm (trial a) 


12.3 nm (trial a) 


4 


64.6 (trial b) 


7.4 (trial b) 


5 


16.9 


9.4 


7 


52.0 . 


7.9 


8 


73.3 


8.2 




Rg. 10: The slicing technique used to accommodate an arbitrary grating 
profile. 



ing the zerotfa diffiraction order. However, for gratings having 
a 1 fim period, the first diffracted order can also be used for 
characterization (or monitoring). Figure 9 shows the results 
of using first order data to predict linewidths of wafer 5. 
Although the PLS results differ from the SEM results at 
higher linewidths, the MMS results are consistent with the 
SEM results for the entire range. The excellent agreement 
between the first order data and experimental results bolsters 
our confidence in the use of die theoretical model to predict 
grating shape parameters, 

VI. COMPLEX STRUCTURE ANALYSIS 

In addition to the primary parameters of interest associ- 
ated with gratings, such as period and linewidth, it is pos- 
sible that more complex structural shapes, such as rounded 
comers and nonrectangular profiles, may be detected and 
characterized by 2-8 scatterometry. Such irregular shapes 



are common iri the production of - submicron structures of 
photoresist since even slight changes in focus and/or expo- 
sure conditions can have a significant effect on the photore- 
sist profile once it is developed. 

We adapted our model to accommodate an arbitrary grat- 
ing profile and generated a Ubrary of scatter signatures from 
complex profiles. The data were generated by slicing each 
profile into many layers of thin rectangular slabs, and then 
ranning the theoretical model over all the layers whOe apply- 
ing the appropriate boundary conditions at each interface. In 
this manner a rounded grating profile is approximated by a 
number of thin rectangular gratings stacked one upon the 
other. This technique has been used in the past to predict the 
dimensions of trapezoidal resist profiles with success.^ 

Two common features encountered in submicron photore- 
sist structures are rounded top and bottom comers. Figure 10 
illustrates how this particular type of profile is sliced into 
different layers. Note that any regions with vertical sidewalls 
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Fig. 9. Prediction results using first order data. 
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Fig. 11. Irregular vs square profiles used in the prediction algorithms^ shown 
here proportionally lo scale. 



need bhly to be represented by one layer. The shape of the 
round portion of each corner is simply a quarter circle' of a 
given radius: A range of corner radii was used iri conjunction 
with the previously predicted nominal linewidth to establish 
the theoretical library. Then the MMS istatistical algorithm 
was implemented to determine the best match between 
theory and experiment, and thus predict the radius of either 
corner; Due to the nature of the chemically amplified resist 
being used, the dimensions of the radii were very small. 
Nevertheless, analyzing profiles with such small dimensions 
(---10 run) provided a rigorous test for the theoretical model 
and prediction algorithm. 

The prediction process for the profiles was done in (wo 
steps. First, the vertical structure library was used to find the 
linewidth of each profile, as previously illustrated: Theii, ust 
ing this predicted linewidth, the MMS algorithm wias run foi^ 
both the top and bottom comer radii of the profiles. By pro- 
ceeding in this manner, we were able to limit our final pa- 
rameter space to just tj;e two comer radii. This drastically 
reduced die computing time required to generate the theoreti- 
cal parameter space. For example, to search over three pa- 
ramelers with ten incremental values per parameter would 
require 1000 theoretical data files. But to limit the search to 
two parameters results in only 100 files. 

Cross-sectional SEMs were perfornied on four of die die 
locations. Three of the locations were used to determine ac- 
tual dimensions related to the overall profile. The fourth Ich 
cation, for which the stepper was defocused intentionally, 
possessed a very complicated profile and was not parameter- 
ized in the usual manner. Instead, the dimensions of the pro- 
file were taken firom the original SEM picture [Rg. 13(a)] 
and implemented into the theoretical model. A scatter signa- 
ture corresponding to this profile was then generated, and 
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Fig. 12. Theoretical scatter signatures for square and irregular shaped pro- 
files compared with the experimental signature. : 



comparisons were made: with the actual experimental data. 
This technique for modeling complex photoresist gratings 
has also been used successfully by Moharam et al}^ 

For this irregular shape, shown as it was mbdeled theo^ 
retically (Fig. 11) in proportion to a square profile, the mean 
square error using 0th order data frorii 0° to 60° incidence 
angle was 12% lower than for the square profile. Thus a 
statistical algorithm would predibt the irregular shape as be- 
ing the more representative one. So, despite the fact that the 
shape has overhanging sidewalls, the stmcture is very well 
characterized by scatterometry. This feature would go unde- 
tected by other techniques apart from cross- sectional SEM 
and the very latest profilometry techniques. Figure 12 shows 
a plot of the theoretical scatter signatures and the experimen- 
tal data taken from the location known to have the irregular 
profile. The particular angular range plotted iaccounts for 
most of the 12% error difference. Clearly the irregular profile 
shows better agreement with the data, illustrating the ability 
of scatterometry to detect this type of shape. 

Regarding the remaining three locations, the dimensions 
of the radii on the top and bottom comers were determined, 
and comparisons were made to SEMs. While the dimensions 
determined by scatterometry are exdct in a computational 
sense, the cross-sectional SEM photos used for comparison 
purposes were difficult to interpret. Because the curvamres. 
of the corners were so small, these features could only be 
estimated froin the SEM pictures using a reticled inspection 
lupe. Although these estimates are subjective and prone to 
error, the comparison shows reasonably good agreement be- 
tween scatterometry and cross-sectional SEM measurements; 



Table IT; Results of complex profile measuiements. 



Exposure/ . 
LxKation . 


Top comer radii 


Bottom comer radii 


Scatterometry 


SEM estimate 


Scatterometry 


SEM estimate . 


30 mj, location A 


0.00 Aim 


No apparent 


0.00 /im 


No apparem 






radius. 




radius - 


60 mJ. location G 


0.09 /jxn 


-b.O4±0.03 /im 


b.OO pm. 


No apparent . 


90 fhJ. location M 








radius 


0.10 /xm 


. -0.68±0.03 Aim 


0.07 /tm 


-0.03 :t 0.03 fjtm 
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Rg. 13. (a) Th6 "inverted S"-shaped profae. (b) the 90 mJ iocarion, with 
both top and bottom rounded corners. 



The results can be seen in Table II. Figure 13(b) shows the 
SEM picture for the 90 mJ location. 

VIL LIMITS OF APPUCATION 

In today's production lines, 0.5 /rni fabrication technol- 
ogy is already standard for high-6nd conunercial semicon- 
ductor devices. In order to keep pace with the industry, it is 
essentiar that the techniques discussed here will be appli- 
cable to dimensions well below 0.5 ^m. 

The performance of these prediction methods on smaller 
dimensions was evaluated by generating theoretical scatter 
signatures for 0.24 /im pitch/0.12 /Am linewidth photoresist 
on Si structures. The illuminating wavelength was 442 nm, 
corresponding to the output of a He-Cd laser. Using the 
resulting angular diffraction signatures, a prediction method 
similar to the PLS algorithm was simulated. The results 
showed that the standard deviation (la) in the prediction was 
5 nm. In addition, we are cunrently investigating actual 0.24 
and 0.35 /tm stnictdres using 633 am illumination wave- 
length. Early results are encouiraging, and it seems reason-^ 
al^le to assume that scatterometry, while using convenient, 
visible wavelengths, v^dll prove viable even as feature sizes 
decrease below 0.1 /tm several years in the future. 

Vm. REPEATABILITY 

The repeatability of an instrument can be defined as the 
variation of measured values when a sample is measm*ed 
repeatedly under constant conditions with that instrument.* 
Repeatability is commonly expressed in multiples of the 
standard deviation, cr, of the repeated measurements. For 
measuring 0.5 /tm features created with a linewidth control 
of 10%, the required 3a measurement repeatability is one- 
third of the linewidth control, or 16 nm. Accuracy is often 
defined as the agreement of a measured value of the "true 
values'* ^'^^ or to an accepted reference value such as a NIST 
standard. Accurate measurements require a set of standards 
for the instrument in question; unfortunately, submicron stan- 
dards are not available in many cases. However, the ability to 
measure accurately is not always as critical as the ability to 
measure repeatably. For example, an instrument with good 
repeatability and linearity can detect changes in die process 
environment. Often, this is all that is required for processing 
applications. 
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To determine the repeatability of the 2-9 scatteix>meterv: 
three experiments were perfomaed, Thd first was designed to 
test the short-term repeatability of thci instrument. In this 
experiirient 30 consecutive zeroth order 2-6 scans were per- 
formed on the saine day on the same part of the Sample (the 
wafer was neVer reni6ved from its mount). For each scan we 
obtained a prediction of the linewidth (using the PLS and 
MMS rbutinesi 60° and 42^ data sets) and calculated the 
mean and standard deviation for the set of iheaisurements. 

The second experiment, which was designed to determine 
the day-to-day reipeatability of the scatterometer, involved 
seven measurements over seven days on the same wafer lo- 
cation. In this experiment the entire scatterometer was started 
cold each day; the sample was removed and repositioned on 
a daily basis. Factors diat would affect this type of repeat- 
ability are mostly due to operator subjectivity when the 
sample is (positioned because of variations across the grating. 

Finally, the third experiment was conducted after the en- 
tire scatterometer optical system was removed, upgraded and 
realigned. This presumably represents the worst possible 
situation in terms of repeatability. The scans were made on 
the same location as in the first two experiments. 

For each experiment the linewidth of the grating was de- 
termined using the standard four prediction algorithms. The 
average and standard deviation associated with each experir 
ment was then calculated. The results, seen in Table III, show 
that the scatterometer demonstrates very good repeatability, 
i.e., a very low standard deviation. For the short-term study, 
the worst 1 cr is 2.0 nm; at best, it is 0.5 nm (using the 42° 
MMS algorithm). For the day-to-day repeatability the worst 
case is slightly higher at 4.0 nm, as expected. The best day- 
to-day repeatability is 0.8 nm, which once again corresponds 
to the 42° MMS algorithm. 

For purposes here the standiard deviation <r is defined to 
be. . . 

V - , (2) 



Recently a new technique for grating liniewidth nieasure- 
ments using an environmental scanning electron microscope 
(ESEM) was repcwted.^^ The ESEM was used in order to 
eliminate the charging effects associated with insulating ma- 
terials such as photoresist when they are used in a conven- 
tional SEM. These charging effects can have an adverse af- 
fect on the linewidth measurement The ESEM sftidy used a 
sample that was comparable to ours: a 1 .0-/im-thick photo- 
resist layer on Si with a nominally 0.5 jjm linewidth. Hie 
ESEM's static Icr repeatability, taken from 10 measurements, 
was 1.3 nm. The day-to-day repeatability was not cited, but 
was presented in a graphical form. From interpolation of two 
of diese data plots we have determined the 1 cr day-to-day 
repeatability to be 1.6 and 2.7 nm, respectively. Thus the 
scatterometric techniques presented; here have repeatability 
figiiies that are as good as, or better than, the latest SEM 
techniques. 
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Table in. Results of the repeaiability study. (Dimensions arc in nanometers.) ;} . . • . ... 



Algorithm Repeatability type Average linewidth . . lo": 

Short term, same day 541.4 2.0 

PLS6O: . Daytoday . . 548J.. . ' 4.0 

Complete rebuild 567.9 n/a 

. ' Short term, same day 551.1 . , ' 2,8 * 

^^PLS42 ■ Daytoday. ■ ■■ 563.9 \ 3.2. 

Complete rebuild. 567.9 . . ■ . , n/a : . 

Short term, same day . . 559.2 * 0.8 ; - 

MMS60 . Daytoday .... 569.7. : : . . . = • 1.0 . . . 

. * Complete rebuild 571.0 . . \. - n/a . 

. . Short term, same day . ... 564.7 ; " . 0.5. . 

MMS42\ . Daytoday , . 565.0 . . . , . . . 0.8 ..... 

. . Complete rebuild . 567.0 . . . n/a . 



IX. CONCLUSIONS 



We have demonstrated that 2-9 scatteroiiietry is capable 
of accurately determining the linewidth of developed photo- 
resist gratings. In addition/ the same techniques can be ap- 
plied to determine the nadre subtle features of the resist lines 
such as rounded comers or overall shape. Comparisons be- 
tween sciatterbmetric and SEM measurements show good; 
agreement for both liiievi^idth and profile measurements. Fur- 
thennore, the precision of the scatterbmeter has been dem- 
onstrated to be excellent, as is evidenced by its low standard 
deviation. In fact, the repeatability of the scatterodaeter is. 
comparable to the best metrology tools available and is sig- 
nificantly lower than a recent, novel SEM technique: 

Characterizing the line shape of developed resist gratings 
is just thie first step in using scatterometjy as a new metrol- 
ogy method. Indeed, in order for scatterometry to be used in 
a seniicondUctor manufacturing environment there are other 
process variables that must be taken into account. For exr 
ample, if there are underlying films that have thickness varia- 
tions from wafer to wafer, linewidth predictions will be af- 
fect^' if these variations are not included in the theoretical 
calibration model. 

At present we are investigating a large, 25 wafer sariiple 
set to determine how additional process parameters affect 
scatterometer measurements. In this experiment, intentional 
variations have been introduced into several of the process 
steps, resulting in a sample set which mimics real process 
drift. Scatterometer measurements are being made for all the 
varying parameters. We will report these results in ia future 
publication. 

ACKNOWLEDGMENTS 

The j authors thank SRC/SEMATECH and tencor/ 
Prbmetrix for their financial support of this work. Specifi- 
cally, they thank I>r. Jimmy Hosch of SEMATECH for pro- 
viding SEMs and valuable discussion and Fat Troccolo of . 
Intel for providing the sample wafers and SEMs. 

'T. R. Corie, SoHc^State Electron: 35, 391 (1992)^ 
^W. L. Stevenson, Proc. SPIE 921. 352 (1988). 
^M. Davidson, Microelectron. Eng. 13, 523 (1991). 
*R. Ucoff, Semicond. InL 17, 56 (1994). 

^ Tolcdo^row. B. W. Smith, Jl K. kbgers, and M. Vaez-lravani. in 



Integrated Circuit Metrology, Inspection, and Process Contml VIH (SPIE, 
BelUngham, WA, 1994):; 

*J. E. Griffith and D. A. Grigg. J. Appl. Phys. 74, R83 (1993). 
'G- Binnig, C. Gcrber, E. Sttrfl, T. R, Albrecht, and C. F. (Juaic, Europhys, 
Lett. 3, 1281 (1987). . 

*S. Manne, H. J. Butt, S. A; C. Gould, and P. K. Hansma, Appl. Phys: Lett 
56, 1758 (1990). ... 

'G. Meyer and N. M. Amer, Appl. Phys. Lett 58; 2100 (1990): 

'**R EHeiz, C. A. Ramos, and P. K. Hansma, J. Vac. Sci. Technol. B 10; 741 

(1992), , . • • 

"D. Nyyssonen, in Integrated Circuit Metrology, Inspection, and Process. 

Contml Vll (SPBE, Bellingham, WA. 1993). pp. 325-335. . 
'^M. T. Postek and D. C. Joy, Solid State TechnoL 29. 145 (1986). 
*'R: R. Henibey and M. B. Weller, Proc. SPIE 1926, 287 (1993). 
*'*J. R. Lowney; M. T Postek. and A. R Vladar, in Ref: 5. pp. 85-96. 
*^K, M. Mohahan, G. Toro-Lira, and M. P. Davidson, in Ref. 11, i^. 336- 

346. .■; ■.. .. 

T- Postek. A. e; Vladar, S. N. Jones, and W. j: Keery, Proc, SPIE 1926, 

268(1993). ! ; 

'^L. A. Firstein and A. Noz, in Integrated Circuit Metrology, Inspection, 

and Process Control V (SPEE, Bellingham. WA, 1991), pp. 81-88. 
''S. S. H. Naqvi;S. H, Zaidi, S. R. J. Bmeck, and J. R. McNeil, J. Vac. Sci. 

Technol. B 12, 3600 ^1994). 
^'^J. R. McNeil et ai, MicroUthog. World 1. 16 (1992). : 
^S. Wttekoek tfr ai, Proc. SPIE 22 (1985). 

^^R. Pforr et al, in Optical/Laser Microlithography V (SPIE, Bellingham, 
WA.. 1992), pp. 594-608. 

^M. P Lang, M.S., University of New Mexico. 1992. " : ' ]■ 

^L. M. Mibier; K. P Bishop, S. S: H. Naqvi, and J. R. McNeil, J. Vac. ScL 
Tfechnol. B II, 1258 (1993). 

^K. C. Hickman c/ oi, J. Vac. Sci. Technol. B IO, 2259 (1992). , 

^J. Snirtevant, S. Holincs, and R Rabidoux. in Advances in Resist technol- 
ogy and Processing DC (SPIE, Bellingham, WA, 1992), pp. 1 14-124. 

^S. S. H. Naqvi et al:, t Opt Soc. Am. A 11, 2485 (1994). 

"S. S. H. Naqvi, J, R. McNeil, R. H. Krukar, and K. R Bishop, Microli^ 
thog. Worid 2, 5 (1993).. . : 

^K. P. Bishop, S. M. Caspar, L. M. Milner, S. S. H. Naqvi. and J. R. . 
McNeil, Proceedings cf the International Conference on the Application 
and Theory of Periodic Structures (SPffi, Bellingham, WA, 1991), pp! 
64-73'. 

^S. S. H. Naqvi, S. M. Caspar, K. t: Hickman, and S: R. McNeil, Appl. 

Opt. 31, 1377(1992).- .. . , . 

Petit, in Topics in Current Physics (Springer, Berlin, 1980). 
^*T. K. Gaylord and M. G. Moharram, Proc. IEEE 73, 894 (1985). 

H. Krukar et al, Proc. SPIE 1926. 60 (1993).. 
"E. V. Thomas and D, M: Haaland, Anal. Cbem, 62. 1091 (1990). . 
^D. M. Haaland and E. V. Thomas. Anal. Cbem. 60, U 93 (1988), . 
^*M. G. Moharam, T. K- Gaylord, G. T. Sincerbox, H. Weriich, and B, 

Yung, Appl. Opt. 23, 3214 (1984). ^ .. . . . 

1986 Annual Book ofASTM Standards (American Society for Testing and 

Materials, Philadclirfiia, 1986, Designation E). 
"T. Yamaguchi. S. Kawata. S. Suzuki. T. Sato, and Y Sato, Jpn. J. Appl. 

Phys. Part 1 32, 6277 (1993). . 



JVST B - Mlcfx>electronics and Nanometer Structures 



